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This article presents an elementary introduction on various aspects of the prototypical integrable 
model the Lieb-Liniger Bose gas ranging from the cooperative to the collective features of many-body 
phenomena. In 1963 Lieb and Liniger hrst solved this quantum field theory many-body problem 
using the Bethe’s hypothesis, i.e. a particular form of wave function introduced by Bethe in solving 
the one-dimensional Heisenberg model in 1931. Despite the Lieb-Liniger model is arguably the 
simplest exactly solvable model, it exhibits rich quantum many-body physics in terms of the aspects 
of mathematical integrability and physical universality. Moreover, the Yang-Yang grand canonical 
ensemble description for the model provides us with a deep understanding of quantum statistics, 
thermodynamics and quantum critical phenomena at the many-body physics level. Recently, such 
fundamental physics of this exactly solved model has been attracting growing interest in experiments. 

Since 2004, there have been more than 20 experimental papers that report novel observations of 
different physical aspects of the Lieb-Liniger model in the lab. So far the observed results to date 
are seen to be in excellent agreement with results obtained using the analysis of this simplest exactly 
solved model. Those experimental observations reveal the unique beauty of integrability. 

PACS numbers: 03.75.Ss,71.10.Pm,02.30.Ik 


INTRODUCTION 

Mathematical principles play significant roles in under¬ 
standing quantum physics. The concept “exact integra¬ 
bility” originally came from the early study of the clas¬ 
sical dynamical systems which were described by some 
differential equations. Usually, the solutions of those dif¬ 
ferential equations are tantamount to the determination 
of enough integration constants, i.e., the integrals of mo¬ 
tions. In this sense, classical integrability is synonymous 
with exact solution. However, conceptual understand¬ 
ing of quantum integrability should trace back to Hans 
Bethe’s seminal work to obtain the energy eigenstates of 
the one-dimensional Heisenberg spin chain with the near¬ 
est interaction in 1931. [l| He proposed a special form of 
the wavefunction — superposition of all possible permu¬ 
tations of plane waves in a ring of size L, namely, 

^ ^ - 

V 

where N is the number of down spins and Vi,..., Vn 
stand for a permutation V of 1, 2, ..., The N\ 
plane waves are A^-fold products of individual exponen¬ 
tial phase factors . Here the N distinct wave num¬ 
bers ki are permuted among the N distinct coordinates 
Xj. Each of the A^! plane waves has an amplitude coef¬ 
ficient which can be in turn determined by solving the 
eigenvalue problem of the Heisenberg Hamiltonian. 

The Bethe’s ansatz appeared to be escaped from physi¬ 
cists’ attention that time. It was only 30 years later 
that in 1963 Lieb and Liniger Q first solved the one¬ 
dimensional (ID) many-body problem of delta-function 
interacting bosons by Bethe’s hypothesis. The exact solu¬ 


tion for the delta-function interacting Bose gas was given 
in terms of the wave numbers ki (i = satis¬ 

fying a set of Bethe ansatz equations, called the Lieb- 
Liniger equations. The spectrum was given by summing 
up all kf. In fact, the Bethe ansatz equations describe 
the roles of individual particles in many-body corrections. 
The significance of quantum integrability is such that the 
roles of individual particles enable one to precisely access 
full aspects of many-body physics for a particular uni¬ 
versity class of many-body systems, for example, spin 
chains, interacting quantum gases, strongly correlated 
electronic systems, Kondo impurity problems, Gaudin 
magnets, etc. Once concerning ensemble statistics, we 
have to properly distinguish physical origin of distin¬ 
guishable and indistinguishable particles. At room tem¬ 
perature, the molecules in the air can be treated as billard 
balls that occasionally collide with each other. The size 
of those molecules is much smaller than the mean dis¬ 
tance between them. The particles are distinguishable. 
However, according to de Broglie matter wave theory, 
the thermal wavelength of a moving particle is given by 
a simple formula: AdB = j(mk^T), where m is 

the mass of the particle, h is the Plank constant, fee is 
the Boltzmann constant, and T is the temperature. The 
temperature becomes very low, the thermal wave length 
increases, and the wave packets start to overlap. Thus, 
the particles are indistinguishable below the degenerate 
temperature. The quantum statistics play an essential 
role under such a degenerate temperature. Below the 
degenerate temperature, there are thus fundamental dif¬ 
ferences between the properties of fermions (with spins of 
1/2, 3/2, ... ) and bosons (with spins of 0, 1, 2, ...). All 
fermions must obey the Pauli exclusion principle, which 
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means that they cannot occupy the same quantum state. 
However, as bosons are not subject to the same restric¬ 
tions, they can collapse under suitable conditions into the 
same quantum ground state, known as a Bose-Einstein 
condensate. The Lieb-Linger model provides an ideal 
ensemble to understand the physics resulting from the 
distinguishable and indistinguishable nature of classical 
and quantum particles. 

Towards a deeper understanding of the physics of the 
Lieb-Liniger gas, a significant next step was made by 
Yang C N and Yang C P on the thermodynamics of this 
many-body problem in 1969.0 They for the first time 
present a grand canonical description of the model in 
equilibrium. In fact, there are many microscopic states 
for an equilibrium state of the system at finite tempera¬ 
tures. The minimization of Gibbs free energy gives rise 
to the so-called Yang-Yang equation that determines the 
true physical state in an analytical way. In the grand 
canonical ensemble, the total number of particles can be 
associated with chemical potential The temperature 
is associated with the entropy S, counting thermal dis¬ 
order. This canonical Yang-Yang approach marks a sig¬ 
nificant step to the exact solutions of finite temperature 
many-body physics. It shows the subtlety of vacuum 
fluctuation, interaction effect, excitation modes, critical¬ 
ity, quantum statistics, thermalization, dynamics, corre¬ 
lations, and Luttinger liquid. 

The name Bethe’s hypothesis was coined by Yang C 
N and Yang C P in the study of the Heisenberg spin 
chain. The Bethe ansatz is now well accepted as a syn¬ 
onym of quantum integrability. ^From solving eigenvalue 
problem of spin-1/2 delta-function interacting Fermi gas, 
Yang0 found that the many-body scattering matrix can 
be reduced to a product of many two-body scattering 
matrices, i.e., a necessary condition for solvability of the 
ID many-body systems. The two-body scattering ma¬ 
trix satisfies a certain intertwined relation, called Yang- 
Baxter equation. This seminal work has inspired a great 
deal of developments in physics and mathematics. The 
Yang-Baxter relation was independently shown by Bax¬ 
ter as the conditions for commuting transfer matrices in 
two-dimensional statistical mechanics. 0, 0 For such ex¬ 
actly solved models, the energy eigen-spectrum of the 
model Hamiltonian can be obtained exactly in terms of 
the Bethe ansatz equations, from which physical prop¬ 
erties can be derived via mathematical analysis. The 
Lieb-Liniger Bose gas0 and Yang-Gaudin modeljl, 0] 
are both notable Bethe ansatz integrable models. 

Yang-Baxter solvable models have flourished into ma¬ 
jority in physics since last 70’s. Later it turned out that 
the Yang-Baxter integrability plays an imp ortant role 
in physical and mathematical studies. |8l4l5l| The Bethe 
ansatz approach has also found success in the realm of 
condensed matter physics, such as Hondo impurity prob¬ 
lems, BCS pairing models, strongly correlated electron 
systems and spin ladders, cold atoms, quantum optics. 


quantum statistical mechanics, etc. The Yang-Baxter 
equation has led to significant developments in math¬ 
ematics, such as 2D conformal field theory, quantum 
croups, knot theory, 2D statistical problems, lattice loop 
models, random walks, etc. Recent research showed that 
there exists a remarkable connection between conformal 
field theory and Yang-Baxter integrability of 2D lattice 
models. Remarkably, conformal field theory has led to 
the theory of vertex operator algebras, modular tensor 
categories, and algebraic topology in connection to new 
states of matter with topological order, such as frac¬ 
tional quantum Hall effect, topological insulators, etc. 
In this elementary introduction to the exactly solvable 
Lieb-Liniger model, we will discuss the rigorousness of 
mathematical integrability and the novelty of quantum 
many-body effects that comprise the beautiful cold world 
of many-body systems in the lab. The content of this ar¬ 
ticle involves understanding the fundamental many-body 
physics through the model of Lieb-Liniger Bose gas. 

The paper is organized as follows. In Section 2, we 
present a rigorous derivation of the Bethe ansatz for the 
Lieb-Liniger Bose gas. In this section we show how a 
field theory problem reduces to a quantum-mechanical 
many-body systems following the Lieb-Liniger’s deriva¬ 
tion. In Section 3, the properties of the ground state 
are discussed, including ground-state energy, excitations, 
cooperative and collective features, Luttinger parameter, 
etc. In Section 4, we introduce the Yang-Yang grand 
canonical approach to the finite temperature physics of 
the Lieb-Linger gas. In this section, we demonstrate how 
the Yang-Yang equation encode the subtle Bose-Einstein 
statistics, Fermi-Dirac statistics, and Boltzmann statis¬ 
tics. In particular, we give an insightful understanding of 
quantum criticality. In Section 5, we briefly review some 
of recent experimental measurements related to the Lieb- 
Liniger Bose gas from which one can conceive the beauty 
of the integrability. 

BETHE ANSATZ FOR THE LIEB-LINIGER 
BOSE GAS 

Wavefunction 

We start with introduction of canonical quantum Bose 
fields satisfying the following commutation rela¬ 

tions: 

[ip{x),^Hy)] = S{x-y), 

[ip{x),tp{y)] = = 0 . 

The Hamiltonian of the ID single component bosonic 
quantum gas of N particles in a ID box with length L is 
given by[2| 

^‘2 p L pL 

H = — dx ^ / dx P, (1) 
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where m is the mass of the bosons, giD is the coupling 
constant which is determined by the ID scattering length 
giu = —2h'^/maiu- The scatteri ng leng th is given by 
aiD = (—aj_/2a„) [1 — C (a^/a i Ij. flbl - flq Here the nu¬ 
merical constant C ~ 1.4603. The model m presents 
the second-quantized form of the Lieb-Liniger Bose gas 
with contact interaction. Q 

In order to process Lieb and Liniger’s solution, we first 
define the vacuum state in the Fock space as tp{x)\0) = 
0, a: G K with (0 | 0) = 1. The equation of motion 
for the field is given by the Heisenberg equation 

idti^{x) = [H,il:{x)]. It follows that the corresponding 
equation of the motion for this model reads 

idt'ipix) = —d‘l'ij){x) + 2c'tp\x)i^{x)'il){x). (2) 

Considering ijj{x) as a classical field, this equation of mo¬ 
tion reduces to a non-linear Schrddinger equation of the 
classical field theory. Moreover, it is easy to show that 
the particle number operator N and the momentum op¬ 
erator P 

^ = j ^ = |['5a;,^'''(a;)]'iA(a;)|dcc(3) 

are commutative with the Hamiltonian ([IJ, i.e., [H, IV] = 
0 and [H,P] = 0. They are among the conserved quan¬ 
tities of this model. The eigenfunction of the A^-particle 
state \P) for the operators iV, iV, and P is given by 

m = P{x)\x), (4) 

where x = {xi,X 2 , ■ ■ ■, a:: at} and 


the dilute quantum gases, when the average distance be¬ 
tween particles is much larger than the scattering length, 
the s-wave scattering between two particles at x^ and 
have the following short distance behavior: [l^ 

r{0+) - r{0-) = -^[‘^'(0+) + tf'(O-)], (7) 

where W(x) is the relative wavefunction of the two parti¬ 
cles and X is the relative distance between the two parti¬ 
cles, i.e., X = Xrj — x^. In the above equation, the prime 
denotes the derivative with respect to x. 

In the model ([6|), the interaction only occurs when 
two particles contact with each other. Following the 
Bethe ansatz,[lj we can divide the wavefunction into 
A^! domains according to the positions of the particles 
0{Q) : XQj < xqj < ••• < xqj,, where Q is the per¬ 
mutation of number set {1,2,..., IV}. The wavefunction 
can be written as ^{x) = '^q 0{Q)4’q{x). Considering 
the symmetry of bosonic statistics, all the ip in different 
domain Q should be the same, i.e., ipQ = ipi, where we 
denote the unitary element of the permutation group as 
1 = (1,2,..., IV}. Lieb and Linger wrote the wavefunc¬ 
tion for the model ([6|) as the superposition of IV! plane 
waves 0 

ipi^Yl ( 8 ) 

V 

where k's are the pseudo-momenta carried by the parti¬ 
cles under a periodic boundary condition. 

Indeed, after solving Schrddinger equation HP{x) = 
EP{x), we can obtain the same s-wave scattering bound¬ 
ary condition © that provides the two-body scattering 
relation among the coefficients A{P) 


\x) = P{xi)ip\x 2 ) ■ ■ ■ iP\xn)\ 0 ). 

Here xj is the coordinate position of the }-th particle. For 
the bosons, the first quantized wavefunction W is symmet¬ 
ric with respect to any exchange of two particles in space 
X, namely, 

!F(...,X5,...,x^,...) = tF(...,x^,...,x^,...). (5) 


A{P') 


kvj 


k-Pj+i 

kvj+i 


+ ic 


— ic 


AiP), 


(9) 


where P' is a permutation obtained by exchanging Pj and 
Pj+i, i.e., P' = {Pi,... ,Pj-i,Pj+i,Pj,Pj+ 2 , ■ • ■ ,Pn}- 
The scattering process in Eq. implies that any two- 
body scattering with momenta fca and kh leads to an anti¬ 
symmetric phase shift A{P') = 


In the following discussion, we set h = 2m = I 
and c = mgiY)/h^. After some algebra, one can find 
that the eigenvalue problem of the Schrddinger equa¬ 
tion iL|'F) = E\P) in IV-particle sector reduces to the 
quantum-mechanical many-body problem which is de¬ 
scribed by the Schrddinger equation E['E{x) = EE{x) 
with the first quantized form of the Hamiltonian 


N 


2=1 


dxj 


N 

2c S{xi — Xj). 

i<j 


( 6 ) 


This many-body Hamiltonian describes N bosons with 
(5-function interaction in one dimension, called the Lieb- 
Liniger model. This is a physical realistic model in quan¬ 
tum degenerate gases with s-wave scattering potential. In 


0(fca 


kh) = 2arctan 



( 10 ) 


When c ^ 0, all the quasi-momenta are different. If there 
are two equal pseudo momenta fca = fcb, we can prove 
that the wavefunction '0i = 0. 

In general, equation ([5]) gives the two-body scattering 
matrix S, A{P') = SA(P) for a quantum many-body 
system. YangQ proved that, if the two-body scattering 
matrix S satisfies the following equation: 

^12(A - m)-5i3(A)43(m) = 523(m)-5i3(A)^12(A - mI}1) 


then the system is integrable. This relation was indepen¬ 
dently found by Baxter(^,y in studying two-dimensional 
statistical models. Nowadays, equation (HU is called 
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Yang-Baxter equation. The Yang-Baxter equation guar¬ 
antees that the multi-body scattering process can be fac¬ 
torized as the product of many two-body scattering pro¬ 
cesses. This factorization reveals the nature of the in- 
tegrability, i.e., no diffraction in outgoing waves. For 
the Lieb-Linger model S matrix is a scalar func¬ 
tion so that the scattering matrix satisfies the Yang- 
Baxter equation trivially. In the scattering process from 
1 to V, the multi-body scattering matrix is defined by 
A{V) = S{Vk)A{l). By using Eq. ([9]), the multi-body 
scattering matrix of this model is given by 


sivk)= n 

Vj<Vi 


kv, 


k-pi — ic 
k-Pi + ic 


( 12 ) 


With the help of Eq. m. the eigen wavefunction of the 
system is given by 


Q,y 


n 

yj<yi 


kv, kv, 
k-p. - fc-p, +icj 


where {7^} are the quantum numbers of pseudo momenta. 
If N is odd, these quantum numbers are integers, whereas 
they are half odd integers when N is even. Eor a given 
set of quantum numbers {7i}, there is a unique set of real 
values {ki} for c > 0. These quantum numbers are in¬ 
dependent of coupling constant c. The total momentum 
can be expressed as 

2 N 

P= — =277r/L, £ = 0,±1, ±2, ... (17) 

i=l 

Eor the ground state, P = 0, where all quasi-momenta 
{fci} are located in an interval {—Q,Q)- Here Q is the 
cut-off. Eor the ground state, the quantum numbers are 
given by 

j = l,...,N. (18) 

GROUND-STATE ENERGY, EXCITATIONS, 
AND CORRELATIONS 


Bethe ansatz equations 


Submitting the periodic boundary conditions 
= 0 ,...) = e'“<7'(..., = L,...) into the 

wavefunction (USD, we can find that the pseudo momenta 
ki satisfies the following Bethe ansatz equations (BAE): 





ki — kj + ic 
ki — kj — ic ’ 


f = l,2,...,iV,(14) 


In the thermodynamic limit, i.e., N, L ^ oo and the 
particle density n = N/L is a constant, the BAE (fTCl) 
can be written in the integral form 

1 

P{k) = — + a{k-q)p{q)dq, \k\<Q, (19) 

J —Q 

where p{k) is the density distribution function of the 
quasi-momenta defined by the particle numbers in a small 
interval of (fc, k -|- Afc), i.e.. 


which are called the Lieb-Liniger equations. When a = 
TT, the wavefunction is anti-periodic; while when a = 0 , it 
is periodic. In the following discussion, we only consider 
the periodic boundary conditions. 

Since partial number N and momentum P are con¬ 
versed quantities of the Lieb-Liniger model, the Hamil¬ 
tonian together with N and P can be simultaneously 
diagonalized. Eor the eigenstate m, the corresponding 
particle number {N) = N. Eor a given set of quasi¬ 
momenta {kj}, the total momentum and the energy of 
the system are obtained as 

N N 

P=iP)=J2 ky E={H)=Y^ kl (15) 

3 3 

The solutions to the BAE (fT4l) provide complete spec¬ 
tra of the Lieb-Linger model. The physical solutions to 
the Bethe ansatz equations require that all the pseudo 
momenta are distinct to each other. The BAE m can 
be written in the form of phase shift function 9{k) as 



p{k) 


1 


lim , , . 
— ^0 T/ /\ k 


Here the cut-off Q is the “Eermi point” of the pseudo 
momenta. It is determined by the particle density n = 
J-Q p{k)dk. Energy per length can also be written as 

^ = [ p{k)k^dk. ( 20 ) 

E J-Q 

For the ground state, the energy depends on the dimen¬ 
sionless scale 7 = Lc/N. Let us make a scaling transfor¬ 
mation 


k = Qx, c = QX, p{Qx) = g{x), (21) 


we find the ground-state energy per particle. 


E 

N 


2m 


eoil), 


( 22 ) 


with 60 ( 7 ) = ^ J^^x'^go(x)dx. Here the distribution 
function go{x) is determined by 


9o{x) 


1 

27T 



9o{y) 

+ {x - uY 


dy. 


(23) 
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with the cut-off condition 7 go(x)dx = A. Equation 
(1^51) is the standard inhomogeneous Fredholm equation, 
which is well understood in mathematics. The Fred¬ 
holm equation can be numerically solved. We will further 
study the ground-state energy later. 

Ground state: from cooperative to collective 

For the ground state, the competition of kinetic energy 
and interaction energy is represented by the dimension¬ 
less coupling strength 7 . When 7 = 0, the system is free 
bosons, and all the particles are condensed at the zero 
momentum state; while if there is a small coupling con¬ 
stant, all the fc’s are distinct. At the limit 7 = 00 , the 
strong repulsion makes the quasi-momentum distribution 
the same as the one of the free fermions. 


The solution of this differential equation is a Hermite 
polynomial and qj is the root of the corresponding poly¬ 
nomial of degree N, i.e., = 0. We set the order 

of qj as qj < qj+i, then we have {2N -I- 1 — q^Y^^ > 
■7T/{qj+i — qj) > {2N+l — q'j^Y^^^. The distribution func¬ 
tion, p(k) = lim^ is thus determined 

^—1 ’ + ) 

hyM 



where the cut-off Q = 2^Jnc is obtained by /_q p{k)dk = 
n. We can see that the quasi-momentum distribution 
function satisfies the semi-circle law where the cut-off 
is the radius of this circle (see Fig. j(a)). The leading 
order of energy (|24ll was also found based on the pseudo 
momentum distribution (|27)l . 


Weak coupling limit: semicircle law 


It is very insightful to examine the physics of the model 
in weak coupling limit. In this case the Bethe ansatz 
roots comprise the semi-circle law. Gaudin firstly no¬ 
ticed such a kind of distribution,]^ followed by several 
groups. 2l|, [22| The Fredholm equation (|23)) is known as 
the Love equation for the problem of the circular disk 
condenser. By using Hutson’s method, Gaudin found 
the density distribution function and energy density 

Q . 2 \ 1/2 


50 (i) 


27rc 


{1 - e) 




1 t 


IbTreQA 
In- 

c / 


e = n 


3/2 


Stt 


(24) 


This result coincides with the perturbative calculation by 
using the Bogoliubov method. 

This problem can also be solved from the original BAE 
(fTTl) . In the weak coupling limit, i.e., Lc <C 1, the quasi¬ 
momenta kj are proportional to the square root of c and 
c/{kj — ki) is a small value. Up to the second order of c, 
the BAE (ITTl) is expanded asjl^ 


1 


N 


(25) 


where qj = kj^jL/2c. If we define a function HN{q) = 
n^i(<Z — 9 i)) we can find that F{qj) = 0 for the polyno¬ 
mial F{q) = H'Y{q) — 2qHYiq)- Both F{q) and Fl^^q) 
are the polynomial of degree N and F{qj) = F[N{qj) = 0, 
and F{q) and FlN^q) are proportional to each other. 
At the large order of g, H{q) = q’^ + ■ ■ ■ and F{q) = 
—2Nq^ + ••■, so that we have F{q) = 2NF[f^{q). It 
follows that 


H'Yiq) - 2qHY{q) + 2NHN{q) = 0. (26) 



Fig. 1. Densities and dressed energies of pseudo momenta 
for the ground state, (a) Solid lines: the dimensionless den¬ 
sities of the pseudo momenta, p{k) = p{k)/c obtained from 
Eq. ([2311; dotted lines: the corresponding dimensionless hole 
densities, ph(^) = Ph.{k)/c. When the coupling strength is 
small, the distribution function p meets a semi-circle law 
For the strong coupling limit, i.e., 7 1, the distribution 

function gradually becomes flatter and flatter, and approaches 
p{k) l/27r. (b) Dimensionless dressed energy is defined by 
e{k) = £{k)/c^, which is obtained from the dressed energy 
equation 1641 . 


Strong coupling limit: fermionzation 


In the strong repulsion limit, the gas is known as the 
Tonks-Girardeau (TG) gas. In realistic experiment with 
cold atoms, it is practicable to observe the quantum de¬ 
generate gas with the strong coupling regime. 2J, 12511 In 
the regime 7 —>■ oo, the Bose-Fermi mapping method [l9| 
can map out the ground-state properties of the Bose 
gas through the wavefunction of the non-interacting 
fermions. For finitely strong interaction, 7 ^ 1 , we can 
obtain the ground-state pseudo momenta from the Bethe 


yrs 

asnatz equations (ICT) [26l 127 1 
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47r^ 

3c3L4 




where Ij take the ground-state quantum numbers Ij = 
— +j — 1, j = 1,..., N. With the help of Eq. (l28l) . 
the ground-state energy of the strong repulsive gas is 
given by 


E 

TT^ 3 

r. 4 

12 

32 


- ^ 

s —n^ 

1 —h 

— 

4 “ —7 

— 1 

L 

3 

7 

72 

q /3 \ 

Vl 5 JJ\ 


■ (29) 


This asymptotic result fits well with the numerical re¬ 
sult obtained by solving the integral BAE (|TO)) (see 
Fig. j(a)). The ground-state energy (1^ can also be ob¬ 
tained from the integral BAE (IT^ by strong coupling 
expansion method. [2^ [2^ We see that for 7 —>■ oo the 
leading order of the ground-state energy reduces to that 
of the free fermions ef = 7 r^n^/ 3 . One can also calcu¬ 
late other important properties such as compressibility, 
sound velocity, and Luttinger parameter, once we know 
the dimensionless function 60 ( 7 ). 

Very recently, Ristivojevic obtained high-precision 
ground-state distribution function for the strong coupling 
regime [ 2 ^ 


p{k) 


1 1 
27r tt^A 

e / 1 

Q2 \^7r2A3 


2 12 - TT^ 8 - 2tt^ 

tt^A^ Stt'^A^ tt^A^ 


+ p)p)+0(V'). (30) 


The ground-state energy per unit length and particle den¬ 
sity are given by 


E_X^f 2 4 120 - 287r2 

L Sttc^ \ ttA tt^A^ IStt^A^ 

80 - 267r2\ 

^ j’ 

_ A / 2 4 24- 47r2 

CTT \ ttA tt^A^ Stt^A^ 

48- 147r2\ 

^ j’ 


(31) 


(32) 


where A = (7 -|- 2 )/ 7 r — 47 r/( 37 ^) -I- 167 r/( 37 ^). The result 
eo = E/{Ln^) is plotted in Fig. „ 


Fig. 2. The ground state energy and sound velocity for dif¬ 
ferent coupling strength 7 . (a) The ground energy density is 
a monotonic increasing function of 7 when the linear density 
is fixed. Black line: the exact numerical result from Eq. (|23]l. 
Red line: the weakly coupling expansion result |[24ll . Blue line: 
the result obtained by the strongly coupling expansion ED. 
(b) The sound velocity is a monotonic decreasing function of 
7 . In the weakly coupling limit, the ratio turns to be infi¬ 
nite, while in the strong coupling limit, Vs/c -^= (see 

Eq. 1(451). 


Elementary excitations: collective motion 


As discussed in previous section, the eigenstates of the 
model are described by the quantum numbers {Ij}. We 
define Ij as a function of the psudomomenta, i.e., I{kj) = 
Ij , thus we have 


m 

L 


A 

27r 



0{k - kj) 
27T 


(33) 


Usually we define the occupied Bethe ansatz roots as par¬ 
ticles, while the unoccupied roots as holes. For a hole 
quasimomentum fch, the quantum number /h = /(fch) S 
Z-|- (iV + 1)/2. We demonstrate the quantum numbers Ij 
for the ground state and excited states in Figs. j(a)-[{d). 
For the ground state, there is no hole below the quasi- 
Fermi momentum, i.e., |/c| < Q, whereas the quantum 
numbers /’s for the holes site outside the interval [Ii,In] 
(see Fig. j;a)). 

In the thermodynamic limit, the pseudo momentum 
distribution function po{k) for the ground state is deter¬ 
mined by Eq. (TT^ with the cut-off of Q. The ground- 
state energy and the particle density can be regarded as 
the function of the cut-off Q, i.e., e{Q) and n(Q). The 
changes over the configuration of quantum numbers for 
the ground state give rise to excited states. Figure j(d) 
shows such an excitation where the particle with a quasi¬ 
momentum fch below the pseudo Fermi point is excited 
outside the Fermi point with a new quasimomentum kg. 
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Fig. 3. Schematic diagrams of the ground state and elementary 
excitations, (a) Quantum numbers for the ground state. The quan¬ 
tum numbers for the ground state are symmetric around the the 
origin. The largest quasimomentum denotes the “Fermi points” 
±Q. (b) Configuration of adding a particle near the right Fermi 
point with the quantum number /at+i, so that the total number of 
particles is -|- 1. (c) A hole excitation. The hole at 7h is created 
so that the total number of particles is A — 1. In panels (b) and (c), 
the parities of their quantum numbers are changed from half-odd 
(or integers) to integer (or half-odds) due to the changes of particle 
numbers, (d) A single particle-hole excitation. A particle at the 
position /h is excited out of the pseudo Fermi sea. In this case, 
the total number of particles is still TV, and the parity of quantum 
numbers dose not change. 


For a single-particle excitation, we decompose the 
pseudo momentum density into two parts, pt{k) = p{k)-\- 
j^5{k — fee), where the delta function term attributes to 
the excited particle and p{k) stands for the density be¬ 
low the cut-off pseudo momenta, i.e., |fc| < Q. In the 
thermodynamic limit, p{k) satisfies the following integral 
equation: 

1 1 

p{k) = 7 . -f / a{k — k')p{k')dk'+—a{k — ke){?)‘^) 

27!' J-Q Tv 

Here we choose a state which satisfies Eq. (IT^ with cut¬ 
off Q as a reference state, i.e., 

I 1 

Po{k) = - -h / a{k — k')po{k')dk'+—a{k — ke^35) 

2tt J_q L 

The addition of the excited particle leads to a collective 
rearrangement of the distribution of the pseudo momenta 
of the N — 1 particles. We denote the difference of the 
pseudo momentum distribution functions between the ex¬ 
cited and the reference state as Ap{k) = p{k) — po{k), 
where poik) is the pseudo momentum distribution of the 
reference state within |fc| < Q. By comparing the inte¬ 
gral BAE (IMll and (1191) , we can find that Ap satisfies the 
following equation: 

1 

Ap{k) = -ra{k — fee) + / ci{k — k')Ap{k'). (36) 

^ J-Q 

The addition of the particle also leads to a shift of the 
cut-off over the psuado-Fermi point of the ground state 
of Wparticle, i.e., AQ = Q — Qq. Here, Qg is the cut¬ 
off for the ground state with the particle number N = 
Lug^Qq), and Q is that for the excited state with N = 
LnciQ) + L J^gdkAp^k) + 1. In the thermodynamic 
limit, AQ is a small value. By comparing the formulas 
of particle numbers between the ground state and the 
excited state, AQnGiQc) = ~[I-q Ap{k)dk+1/L]. The 
prime denotes the derivative with respect to the the cut¬ 
off Q. The corresponding excited energy 

fQ 

AE(fcc) = L Ap{k)k^dk + kl + Le'{Q)AQ 

J-Q 

rQ 

= L Ap{k){k'^ - p)dk + (kl - p), (37) 

J-Q 


where p is the chemical potential, and p = dE/dN = 
SG{QG)/nG(QG)- Later we will further prove that the 
dressed energy can be expressed as 

rQ 

£(k) = k^ — p+ / a(k — k')£(k'). (38) 

J-Q 

For convenience, we introduce a useful relation between 
the dressed energy and density. Assuming that two func¬ 
tions f{k) and g(k) satisfy the following equations: 


rQ 

f{k) = foik)+ a(k-k')f{k')dk', 
J-Q 


rQ 

9{k) = 9o(k) + a(k- k')g{k')dk', 
J-Q 


\k\ < Q, 

\k\ < Q, (39) 


where fo(k) and go{k) are the driving terms of these in¬ 
tegral equations. Then we have the useful relation 

rQ rQ 

/ fik)go{k)dk = / g{k)fo{k)dk. (40) 
J-Q J-Q 

By using Eq. (l40ll . from Eqs. (IM)) and (1^ . we thus prove 
that the dressed energy is nothing but the excitation en¬ 
ergy of a particle, 

e(fce) = kl - p + J^Q a(ke - fc')e(fc') 

= kl-p + L r (kl - p)Ap(k') = AE(ke). (41) 
J-Q 


On the other hand, the hole excitations impose an ad¬ 
ditional condition, namely, the maximum momentum is 
riTT. Similarly, one can prove that the excited energy 
reads 


AF;(fch) = -e(fch). (42) 


We plot the dressed energies for different coupling 
strengths in Fig. j(b). In order to see clearly excitation 
spectra, we need to calculate the total momentum us¬ 
ing Eq. (ED. For both the particle-hole excitation and 
the Type-H hole excitation, the total particle numbers of 
particles do not change, as shown in Figs. j(d) andj(b). 
Therefore, the total momenta of the excited states are 
given by 

r^e,h 

P = 7111 — 27: / po(k)dk. (43) 

Jo 

For the low-energy excitations, fce — Q is a very small 
value. Therefore, the low-lying behavior can be described 
by a linear dispersion relation 


AE = eG(Q)\k ± Q\ = Vs\P\, 


£g(Q) 

27rpG(Q) ’ 


(44) 


where Vs is the sound velocity of the system. In the above 
equation, ± corresponds to the excitations at left/right 
Fermi point, respectively. Nevertheless, for |fc| > Q, the 
relation (H^ is the dispersion relation for particle-hole 



excitations. This linear spectra uniquely determine the 
universal Luttinger liquid behaviour at the low tempera¬ 
tures. The essential feature resulted from the dispersion 
relation (1441) is such that the system exhibits the confor¬ 
mal invariant in low energy sector. We will further dis¬ 
cuss this universal nature of the ID many-body physics. 
For the strong coupling, we can obtain the sound ve¬ 
locity from the ground-state energy through the relation 

= VS?’ namely, 


Bethe ansatz equations o, near c —> ±oo, the com¬ 
pressibility is given by 


— = 27r^n — 

K 


IGtt^ 2 

- n 

c 



+ 



However, for the Lieb-Liniger gas with weak repulsive 
interaction (0 < c <C 1 ), compressibility is given by 


1 

At 


= 2c-^c3/2. 
Tryn 


(47) 


Vs 


2'Kn 



12 

712 



(45) 


In Figs. d)a) andj^b), we present the analytical and nu¬ 
merical results for the dimensionless energy 6 ( 7 ) = e/v? 
and Vs/c, respectively. 



Fig. 4. Dispersion relations in elementary excitations: adding 
one particle and adding one hole excitation, (a) Solid lines: 
the spectra for adding one particle to the ground state corre¬ 
sponding to the configuration presented in Fig. j^b). Dashed 
lines: the dispersion relations for the hole excitation, see the 
configuration shown in Fig.|{c). For 7 = 1 , the thin blue line 
is calculated by using Eq. Hill . We see that both the parti¬ 
cle and hole excitations comprise the dispersion relations with 
the same velocity, (b) The particle-hole excitation spectra, 
see Fig. ((d). The linear dispersion relation is seen for the 
long-wavelength limit, i.e., the momentum tends to be zero. 


The super Tonks-Girardeau gas-like phase 


In regard of the strong interacting bosons in ID, the 
super Tonks-Girardeau gas is particularly interesting. It 
describes a gas-like phase of the attractive Bose gas which 
was first proposed in a system of attractive hard rods 
by Astrakharchik et aL[30j Batchelor and his coworkers 
showed its existence of such a novel state in the Lieb- 


Liniger model with a strong attraction. ^ Due to the 
large kinetic energy inherited from the repulsive Tonks- 
Girardeau gas, the hard-core behavior of the particles 
with Fermi-like pressure prevents the collapse of the su¬ 
per TG phase after the switch of interactions from repul¬ 
sive to attractive interactions. 32-3^ In fact, the energy 
can be continuous in the limits c —>■ ±oo. ,(,From the 


In contrast, the energy for the gas-like phase (super 
Tonks-Girardeau gas) in weak attractive interaction 
limit (—1 <C c < 0 ) is given by 

E=^n^-\c\n^. (48) 

Thus, the compressibility of the super Tonks-Giraread 
gas with weak attraction is given by 

— = 47 r^n — 2 |c|. (49) 

K 


Such different forms of compressibility reveal an impor¬ 
tant insight into the root patterns of the quasi-momenta. 
We now show the subtlety of the Bethe ansatz roots for 
the super Tonks-Girardeau state [32| in Fig. „ 



k (units of 27 i/L) 


Fig. 5. Quasi-momenta for the ground state of the Lieb- 
Liniger gas with the strong repulsion and the supper Tonks- 
Girardeau gas of the attractive Bose gas for a large value of 7 . 
It is obvious that the super Tons—Girardeau gas has a larger 
kinetic energy than the free Fermion momentum. Thus the 
quantum statistics of the super Tonks-Girardeau gas is more 
exclusive than the free fermion statistics. 


Luttinger parameter and correlation functions 

The ID integrable system gives rise to the power law 
behavior of long distance or long time asymptotics of 
correlation functions for the ground state. The effec¬ 
tive Hamiltonian can be approximately described by the 
conformal Hamiltonian which can be written in terms of 
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the generators of the underlying Virasoro algebra with 
the central charge C = 1. The low-lying excitations 
present the phonon dispersion AE(p) = VgP in the long- 
wavelength limit. In this limit, all particles participate 
in the excitations and form a collective motion of bosons 
which is called the Luttinger liquid. 36| The Lieb-Liniger 


field theory Hamiltonian can be rewritten as an effective 
Hamiltonian in long-wavelength limit, which essentially 
describes the low-energy physics of the Lieb-Linger Bose 
gas 


H = 


da; 


2 2ttK 



(50) 



where the canonical momenta H conjugate to the phase 
(j) obeying the standard Bose commutation relations 
n(y)] = i(5(a; — y). dx4> is proportional to the den¬ 
sity fluctuations. In this effective Hamiltonian, VgjK 
fixes the energy for the change of density. In this ap¬ 
proach, the density variation in space is viewed as a su¬ 
perposition of harmonic waves. 

For example, the leading order of one-particle correla¬ 
tion ('0^(x)'!/'(O)) is uniquely determined by 

the Luttinger parameter K. The Luttinger parameter is 
defined by the ratio of sound velocity to stiffness, namely, 


^3eo(7)-27^ + h^^’ 

where is the sound velocity and wat is the stiffness and 
defined as 


L d^E _ I L2 d2E 
ttH dN'^ ’ ^ V rnN dL'^ 


(52) 


In Eq. (jSTI) . the second expression of the Luttiger pa¬ 
rameter is used for numerical calculation with the help 
of Eq. (HH). 

Using the asymptotic expansion result of the ground- 
state energy ((^ for weak counting and (HH) for the 
strong counting regimes, we find the asymptotic forms 
of the Luttinger parameter K for the two limits 


I ' 

= 7r|^7- ^7^^^ 

4 4 IbTT^ 

^I7>1 — IH-1- 2 - 

7 7^ 6')'^ 


(53) 

(54) 


In Fig. n we show that these asymptotic forms of the 
Luttinger parameters provide a very accurate expression 
throughout the whole parameter space. The correla¬ 
tion functions can be calculated by the conformal field 

I ■ I I I ■ 

theory. [ 


Fig. 6. Luttinger parameter K versus interaction strength 
7 . Blue dashed line: the result obtained from Eq. (154 l l for 
the strong coupling limit; red dashed line: result obtained 
from Eq. (I53I I for weak coupling regime; black solid line: the 
numerical result obtained from Eq. (Eg. The analytical result 
of the Luttinger parameter is in a good agreement with the 
numerical result. 


Usually, the long-distance or long-time asymptotics of 
correlation functions of the ID critical systems can be 
calculated by using the canonical field theory (CFT). [H 
^From the CFT, the two-point correlation function for 
primary fields with the conformal dimensions is given 
by® 

A„-2Tvi{NAD}y/L 

Goiy,T) = ^- - - . ^2^- , (55) 

[VT + lyY^ [VT — lyY^ 


where r is the Euclidean time, v is the velocity of light, 
Goix,t) = (Glb^a;, t)d(0,0)|G') is the correlator for the 
field operators 0(0,0) and 0^{xY)- Equation (1551) in¬ 
volves the contributions of the excited states which are 
characterized by numbers AD, N^, and AN. Here N~ 
(or iV+) describes the elementary particle-hole excita¬ 
tions by moving atoms close to the left (or right) pseudo 
Eermi point outside the Eermi sea with adding N~ (or 
V+) holes below the left (or right) Eermi point. AN is 
the change of particle number and it characterizes the 
elementary excitations by adding (or removing) particles 
over the ground state. N^ and AN are not enough to 
describe all elementary excitations. It is necessary to 
introduce the quantum number 2AD to denote the par¬ 
ticle number difference between the right- and left-going 
particles. ^Erom the analysis of finite-size corrections of 
the BAE, the total momentum and excited energy of the 
low-lying excitations in terms of these quantum numbers 
N^, AN, and D are given by 


O-jT 

AP= — [ANAD + N+ - N-] + 2ADkF, 
L 


AE = 


2ttv 


-{AN/Zf -t- {ADZf + N+ + N- 


(56) 

(57) 


where Z = 2'Kp{Q) is the dressed charge at the pseudo 
Fermi point for the ground state. ^From the conformal 
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field theory, the excited energy and momentum are given 

by 

zTT?^ 

A£;= —w(Z\+ + Z\-), (58) 

1j 

‘lir 

AP= — ^(A+-A-) + 2Ai:ifcF, (59) 


By comparison between the two results obtained from 
finite-size corrections dMl), dSZl), and CFT ([Ml), the 
conformal dimensions are analytically obtained as a func¬ 
tion of AN, AD, and the dressed charge Z as well, 
namely, 

2A=^ = 2N^ ± AN AD + {ADZf + ^{Z-^AN)%m) 


It turns out that only the low-energy excitations de¬ 
termine the long distance or time asymptotic behavior of 
the correlation functions. Given that the operator '(/>(x) 
is a prime operator of this 17(1) symmetric system, the 
correlation functions of the prime field have a universal 
power law decay in distance 


(^^(A2/)^(0,0)) 


exp(2iAi7fcF?/) , . 

{vT + {vT — 


We find that the Luttinger parameter K = Z^ hy 
comparing with the result of the Luttinger theory, 
(^/’'^(j/)'0(O)) ^ For example, in the strong cou¬ 

pling regime, the dressed charge Z = 1 + 2/7 — STr^/Sy^ 
from the dressed charge equation. ll| Submitting it into 
K = Z^, we obtain the same result as Eq. (IMl). 


YANG-YANG THERMODYNAMIGS AND 
QUANTUM GRITIGALITY 


In 1969, Yang C N and Yang C P presented a grand 
canonical ensemble to describe finite-temperature ther¬ 
modynamics for the Lieb-Linger model. [3| The Yang- 
Yang method has led to significant developments in quan¬ 
tum integrable systems. [l2. 1^ U| This approach allows 
one to access full finite temperature physics of the mod¬ 
els in terms of the thermodynamic Bethe ansatz (TBA) 
equations, [l^ In the grand canonical ensemble, we usu¬ 
ally convert the TBA equations in terms of the dimen¬ 
sionless chemical potential fl = y-jc? and the dimension¬ 
less temperature T = Tj(? with the interaction strength 
c. It is also convenient to use the degenerate temperature 
as an energy unit, i.e., Td = j2m. ^ l43| It is very 

insightful to discuss the critical phenomena of the models 
in terms of the dimensionless units. We first discuss the 
Yang-Yang grand canonical ensemble below. 


The Yang—Yang grand canonical ensemble 


tures, any thermal equilibrium state involves many mi¬ 
croscopic eigenstates. As discussed in previous section, 
these eigenstates are characterized by different quantum 
numbers {Ij}, see Eq. (1161) . In the thermodynamic limit, 
I{k) is a monotonic function of pseudo momenta k. 111 
We define dI{k)/{Ldk) = p(k) + ph{k), where ph is the 
density of the holes. dFrom Eq. we can obtain the 
integral BAE for arbitrary eigenstate as 


1 7“ 

p{k) -I- ph{k) = - -1- / a{k — k')p{k')dk'. (62) 

J — CXD 


Here we should notice that integral interval can extend 
to the whole real axis, i.e., particles can occupy any real 
quasimomentum. 

In order to understand the equilibrium states of the 
model, it is essential to introduce the entropy. In 
a small interval dk, the number of total vacancies is 
L[p{k) -(- pii{k)]dk with a number of Lp{k)dk particles 
and a number of Lp\,{k)dk holes. These particles and 
holes give rise to microscopic states 

^ [L{p{k) + pi,{k))dk]l 
[Lp{k)dk ]! [Lph{k)dk ]! 

In the thermodynamic limit, [L{p{k)+phik))dk] ^ 1 and 
dk —^ 0, with the help of Stirling’s formula, the entropy 
in this small interval is given by 

dS" = IndVF « L{pln[l + p] + Ph ln[l -b ri~^]}dk, 

where ri{k) = p\^{k)/p{k). The total entropy is S = f dS”. 
It can be understood from this procedure that the en¬ 
tropy involves the disorder of mixing the particles and 
holes in the pseudo momentum space. For regions with 
zero p(k) or zero Ph{k), no disorder occurs, i.e., dS” = 0. 
Therefore, the entropy for the ground state is zero. It is 
worth noting that the above discussion is valid only in 
the equilibrium state. 

The Gibbs free energy I? is the thermodynamic poten¬ 
tial of grand canonical ensemble for this model. 


n = E -TS - pN, (63) 


where y is the chemical potential, and the particle num¬ 
ber is given by N = L J p(k)dk. In thermal equilibrium, 
the true physical state is determined by the conditions 
of minimizing the Gibbs free energy. Making a virtual 
change 5p, Sp^ in the thermal equivalent state, we take 
the variation of the free energy such that 


Sn = SE- TSS - ySN = 0. 


Here we should notice that the variations Sp and dph are 
not independent in view of the integral BAE (15^ . This 
minimization condition leads to the TBA equations in 
terms of the dressed energy Q 


For the ground state, the set of quantum numbers m 
provides the lowest energy. However, at finite tempera- 


/ CX) 

a(k — k')£-(k')dk', (64) 

-00 
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where the dressed energy is defined by eik) = T\nri{k). 
In the above equation, we also denote e-{k) = —Tln[l + 
Using the TEA equation (1^ . we fur¬ 
ther obtained the grand thermodynamic potential fi = 
^ J e-(k)dk. It follows that the pressure is given by 


P = - 



1 

277 


S-(k)dk. 


(65) 


This serves as the equation of state (EH) from which we 
can calculate the thermodynamics of this model at finite 
temperatures. We can obtain the zero-temperature and 
finite-temperature phase diagrams of the Lieb-Liniger 
model. In the zero-temperature limit, the TEA equa¬ 
tion ([M]) reduces to the dressed energy equation (1551) 
in the limit T —>■ 0. From the standard thermody¬ 
namic relations, one can calculate the particle density 
n = 5^p|c,T) entropy density s = compressibility 

K* = 5^p|c,t, specific heat Cv = Td^p\fj,^c in a straight¬ 
forward way. 



Fig. 7. Quantum critical regimes for the Lieb—Liniger model, 
(a), (b) Dimensionless entropy and specific heat in the T-p. 
plane, respectively. For T 3> 1, it is the HT regime. The TLL 
with the dynamical exponent z = 2 and correlation length 
exponent i/ = 1 lies in the region T <g 1 and fi > fic- For 
T ^ \fi — /Tel, the QC regime with z = 2, d = 1, and u = 1/2 
fans out near the quantum phase transition point pc = 0. It is 
obvious that the entropy and the specific heat have singularity 
properties near the critical point fic = 0. 


In view of the grand canonical ensemble, there exists 
a quantum phase transition at the chemical potential 
/ic = 0 at zero temperature. Universal thermodynam¬ 
ics is expected for the temperature under the quantum 
degenerate regime T < Td, where Td = fi^n^ /2w.[4l. 

At high temperatures (HT), i.e., T ^ Td, the system be¬ 
haves like a classical Eoltzmann gas. For /i < /ic = 0 and 
at low temperatures, the density is very low and the gas 
becomes de-coherent. This phase is semiclassical (SC). 
Whereas for p > ptc and the temperature T < \p — pd, it 
shows the Tomonaga-Luttinger liquid (TLL) phase. The 
quantum critical regime lies between SC and TLL for 
the temperature T :S> \p — pc\. We show such different 
regimes through the entropy and specific heat in Fig. □ 


The Yang—Yang equation and quantum statistics 

Dynamical interaction and thermal fluctuation drive 
Lieb-Liniger model from one phase into another. In par¬ 
ticular, under the degenerate temperature Td, the model 
has three distinct phases: semi-classical, quantum criti¬ 
cal, and the TTL critical phases. At high temperatures, 
there does not exhibit universal behaviour. When the 
temperature tends to be infinity, the system reaches the 
Eoltzmann gas. Therefore, the Yang-Yang equation ((64ll 
encodes different quantum statistics. For example, when 
the coupling strength 7 —>• 0 , the system behaves as the 
free bosons; for the strong coupling limit, it behaves like 
free fermions; at high temperatures, the system becomes 
the Eoltzmann gas. In the following, we rigorously derive 
such quantum statistics in an analytical way. 

When the coupling strength 7 turns to be zero, the 
integral kernel a{k) —>■ 6{x). The dressed energy in this 
limit can be expressed as 

lim£(A:) = Tin [e*-^ -idIT _ (^ 55 ^ 

from which we obtain the thermal potential per length 
limc^o D/L = 2 ^l\Q^^~d)IT _ 1 ] de. Thus, the dis¬ 

tribution function satisfies the Eose-Einstein statistics 

If 7 —>■ 00 , the dressed energy reads 

lim e{k) = k^ — p. ( 68 ) 

c—^oo 

In this limit, the Eethe ansatz equation (1621) naturally 
reduces to the form 

= 27r(I + ph(fc)MA))’ 

which indicates the Fermi-Dirac statistics. Conse¬ 
quently, the thermal potential per unit length is given 
by limc_).oo D/T = + 1] de. This result 

gives rise to the Fermi-Dirac statistics 

We further remark that equations dMll-dlQl) are valid 
for arbitrary temperature. If the system is under the 
quantum degeneracy, the quantum statistical interaction 
is important. Thus, the particles are indistinguishable. 
At high temperatures, the Yang-Yang equation (|64)l gives 
rise to the Maxwell-Eoltzmann statistic such that the 
particles are distinguishable. In the weak coupling limit 
and high-temperature limits, it is very convenient to 
consider Viral expansions with the Yang-Yang equation 
(l64l) . namely, 

Q-dk)/T ^ /T^!T^Aqa(k-q)\nil+Ze->'"^ 
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Here Z = e^^^ is fugacity. After some algebra, we find 
that the pressure up to the second Viral coefficient is 
given by 

7 - 3/2 

P = Po + ^=^^P 2 , (72) 

V ZTT 

where p 2 = dq'a(2q')e~^'^ reveals the two- 

body interaction effect. In the above equation, po = 
— ^ dfcln — Ze~^ is the pressure of the free 

bosons. The result GH) gives the Maxwell-Boltzmann 
statistics in the limit of T —>■ 00 . 

It is remarkable to discover the universal low tem¬ 
perature behavior of the Lieb-Liniger model with the 
Bose-Einstein statistic and Fermi-Dirac statistic. The 
Yang-Yang equation (Ifi4l) provides full physics of the 
Lieb-Liniger model which _goes beyond that can be found 
by Bose-Fermi mapping. [l9| In fact, for strong coupling 
limit, i.e., 7 1 , the system can be viewed as an ideal 

gas with the fractional statistics, When the coupling 
strength is very weak, the ground state behaves like a 
quasi BEC. The Bogoliubov approach is valid in the weak 
coupling limit T <C ^7 l.jS, 0, The result (17^ 
is also a good approximation for the weak coupling Lieb- 
Liniger gas. 


Luttinger liquid and quantum criticality 

Equation of state 


In low-energy physics, T Jdj low-lying excitations 
form a collective motion of bosons. The linear relativistic 
dispersion near the Fermi points results in the TLL be¬ 
havior. At finite temperatures, the TLL can be sustained 
in a region of T < |/i —pd in the T-/r plane (see Fig.p). In 
the TLL phase, we can take the Sommerfeld expansion 
with the TBA equation (IMl) . By iterations, the pressure 
with the leading order temperature correction is given by 


P = Po + 


Po(Q) 

3 e'iQ) 


Po 


ttT^ 
6vs ’ 


(73) 


which gives the free energy per unit length as the field 
theory prediction 


F{T)/L = fin— p 


nCjkBT) 
6hv 


(74) 


with the central charge (7=1. Here we took fce = 1, so 
that in TLL phase the specific heat is linear temperature- 
dependent 


_ ttT 
3us 


(75) 


This is a universal signature of the TLL. 

However, for the temperature beyond the crossover 
temperature, i.e., T > T* ^ \p ~ Pc\, the excitations 


give a non-relativistic dispersion, i.e., AA ~ p^. The 
crossover temperature T* can be also determined by 
the breakdown of linear temperature-dependent relation 
given by Eq. The crossover is also evidenced by 

the correlation length. 45 | 

For strong coupling and low temperatures, i.e., 7 1 

and r <C 1 , the pressure is given by [ill 


r3/2 / A/T\ 


1 + 


7-3/2 

2i/7rc' 


Li3/2( 


(76) 


where A = fi+2pjc+{T^^^j 2 ^c^)Li 3 / 2 {—e^^'^) and the 
polylogarithm function is given by Li„(x) = pr- In 

the above equation, the pressure p gives a close form of 
the equation of state. Using the standard thermodynam¬ 
ical relations, we can analytically calculate the particle 
density, compressibility, and the specific heat 


n = — 


2 v^ 




7-1/2 j 

rp3/2 


1/2 


1 - 


V^( 


7-1/2 j 


1/2 + 


-fl/2 + 2 / 3/2 


(77) 




^T-1/V-1/2 + ^/-l/2/l/2 

f-l/2fl/2 -5'7"/-1/2/3/2 


Cy 

T 


7^3/2j,2 - , 

dT 


(78) 


,A 


2v7r 


/J,,C 


7-2 


80F 20F 


/-1/2 + o 


(79) 


respectively. Here /„ = Li„(—e-^/^). 


Quantum criticality 

At zero temperature, the quantum phase transition 
from the vacuum phase into the TLL at the critical point 
fic = 0 occurs in the Lieb-Liniger Bose gas. Accord¬ 
ing to the renormalized group theory, universal scaling 
properties are expected in the critical regime at low tem¬ 
peratures (see Fig. 0 ). In 2011, Guan and Batchelor in¬ 
vestigated the quantum criticality of the Bose gas and 
found that the equation of state Gi) reveals the uni¬ 
versal scaling behavior of quantum criticality in terms 
of the polylogarithm functions. [ 2 ^ It is straightforward 
from the equation of the state to derive the universal 
scaling form of the density as 

+ (80) 

where the background density uq = 0. The scaling func¬ 
tion iF{x) = — j^Lii/ 2 (~e^) read off the dynamic criti¬ 
cal exponent z = 2 , and the correlation length exponent 
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1 / = 1/2. It is particularly interesting that the finite 
temperature density profiles of the ID trapped gas can 
map out the quantum criticality with these universal ex¬ 
ponents for the Lieb-Liniger gas. The density curves at 
different temperatures intersect at the critical point, see 
Fig. Jl^a). The universal scaling behavior of compressibil¬ 
ity K* is given by 

= > (81) 

where /tg = 0 and Q(x) = — ^;^Li_i/ 2 (a;). This scaling 
function again reads out the dynamic critical exponent 
z = 2, and the correlation length exponent i/ = 1/2. The 
intersection at the critical point for different tempera¬ 
tures attributes to the universal scaling form m (see 
Fig.j:b)). 



A A 


Fig. 8. Universal scaling behaviors of the density and com¬ 
pressibility at quantum criticality, (a) Density shows the 
universal scaling behavior given by Il80|l . (b) Compressibil¬ 

ity presents the universal scaling behavior (EJ. The inset in 
(b) shows the collapse of temperature-rescaled compressibility 
with respect to the argument (A — fIc)/T. 



Fig. 9. Quantum critical behaviour of the entropy and spe¬ 
cific heat, (a) The entropy divided by temperature s = s/c 
has a universal scaling bahavior at the quantum criticality, (b) 
Specific heat divided by temperature Cv = Cv/T has a univer¬ 
sal scaling bahavior at the quantum criticality. The inset in 
(b) shows the collapse of Cv\/t with respect to the argument 
(A-Ac)/^. 


Near the critical point, the specific heat divided by the 
temperature Cy = c^/T obeys the following scaling form: 


rp-l/2 


^ Li3/2( - ^ 






2^T 


Lii/2( - 





2 

Li-1/2 ( 



(82) 


The specific heat at different temperatures has two round 
peaks near the critical point fic = 0. These peaks mark 
the crossover temperatures that distinguish the TLL 
and semi-classical gas phases from the quantum critical 
regime. This is a very robust signature for the existence 
of the crossover temperatures in the ID Bose gas. This 
scaling law of the entropy and specific heat is shown in 
Fig.„ 


The local pair correlations and Tan’s contact 


In the study of the interacting Bose gas in one di¬ 
mension, an important property is the local two-body 
correlation function g 2 . Physically speaking, this func¬ 
tion describes the rates of inelastic collision between pairs 
of particles. This quantity reflects the probability 

that the two particles site at the same size. It is also 
known as contact that strikingly captures the universal¬ 
ity of ultracold atoms. This has been described by Tan’s 
relations. Tan’s contact, which measures the two- 

body correlations at short distances in dilute systems, is 
a central quantity to ultra-cold atoms. It builds up uni¬ 
versal relations among thermodynamic quantities such as 
the large momentum tail, energy, and dynamic structure 
factor, through the renowned Tan’s relations, see recent 
developments in this research. [H^ 

Knowing that the two-body correlation can lead to the 
classification of physically distinct regimes, for example, 
the Tonks-Girardeau regime with 32 —>■ 0, the Gross- 
Pitaevskii regime with 32 = 1 and the very weak coupling 
or fully decoherent regime with g 2 = 2. The two-body 
correlation function can be calculated for these differ¬ 
ent regimes by expression g 2 = (!F'l(a;)>Fl'(a:)lF(x)lF(x)). 
where W is the field operator in second quantization. At 
T = 0, we have dAg/dc = Lg 2 , where Eq is the ground- 
state energy. For weak coupling limit, the local pair cor¬ 
relation 32 /'^^ = 1 — 2 ^ 7 / 71 . For strong coupling limit, 
32 /n^ = 477 ^/ 87 ^ (I — 6 / 7 ). In general, the local pair cor¬ 
relation can be used to study phase coherence behavior 
at finite temperatures. Introducing the free energy per 
particle, i.e., /(y, T) = F/N, the normalized two-particle 
local correlation is defined as 

{E^ix)EHx)Eix)Eix)) 2m fdf{j,T)\^ 

92 = - 7 ^- = 7777 1 - 77 - ) 


/j2^2 Qj 

For strong coupling regime, the local pair correlation 
function is obtained from the equation of state TT], l48l . 

MUl 

^ (84, 

d 


37 7'^47r2r|j 


On the other hand, in one dimension the fundamen¬ 
tal thermodynamic relation in a harmonic trap is given 
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byliil 

dp = nd^ + sdT — — CdaiD, (85) 


where ps and C are the densities of superfluid and contact, 
respectively. In this relation, w = Vs — v^ is the difference 
between the velocity of the superfluid and normal com¬ 
ponents. Maxwell relations build the general connections 
between the contact and other physical quantities such 
as 


[dp 


T,air> 


V^OiD / ^ 2" 

(^) 

ydaiB J 


Furthermore, we can obtain the contact through the fol¬ 
lowing relations: 




2tt 


fl/2f3/2 


X 



( 86 ) 


The Tan’s contact of the Lieb-Liniger gas does not have 
the usual scaling behavior which was found for the in¬ 
teracting Fermi gas in Ref. This is mainly because 
the critical field pc = 0 which does not depend on the 
scattering length oid (see Fig.,). 



Fig. 10. Contact divided by versus chemical potential. 
The critical scaling behavior of the contact for the Lieb- 
Liniger gas is different from that for the ID interacting Fermi 
gas. l28ll 


EXPERIMENTAL DEVELOPMENT RELATED 
TO THE LIEB-LINIGER GAS 


Over the past few decades, experimental achievements 
in trapping and cooling ultra-cold atomic gases have re¬ 
vealed beautiful physics of the cold quantum world. In 


particular, recent breakthrough experiments on trapped 
ultracold bosonic and fermionic atoms confined to one 
dimension have provided a precise understanding of sig¬ 
nificant quantum statistical and strong correlation ef¬ 
fects in quantum many-body systems. The particles in 
the waveguides are tightly confined in two transverse 
directions and weakly confined in the axial direction. 
The transverse excitations are fully suppressed by the 
tight confinements. Thus, the atoms in these waveg¬ 
uides can be effectively characterised by a quasi-lD sys¬ 
tem. Thus, ID effective interaction potentials can be con¬ 
trolled in the whole interacting regime by the underlying 
3D scattering with tight confinements in the two trans¬ 
verse directions. |16l4l8l | In such a way, these ID many- 
body systems ultimately relate to the integrable mod¬ 
els of interacting bosons and fermions. It is now possi¬ 
ble to realize effectively one-dimensional quantum Bose 
gases, in which the interaction strength between ultracold 
atoms is tunable, see recent reviews, [s^, El These ex¬ 
periments have successfully demonstrated the anisotropic 
confinements of atoms to one dimension by optical waveg¬ 
uides, see a feature review article. Particularly strik¬ 
ing examples involve the measurements of momentum 
distribution profiles. 1^ the ground state of the 
Tonks-Girardeau gas,[^ quantum correlations. 50 - 161 
Yang-Yang thermodynamics, 62|, l63| the super Tonks- 
Girardeau gas, 641 quantum phonon fluctuatio ns. 651 467 
elementary excitations and dark solitons, 68], l69l| ther- 
molization and quantum dynamics. 7^721 More experi¬ 
mental developments of the Lieb-Liniger model are listed 
in Table 


Table 1. Experiments of Lieb-Liniger gas. 


quantum dynamics 
thermalization 
solitons 
fermionization 
YY thermodynamics 
strong coupling 
phase diagram 
3-body correlations 
excited state 


®’'Rb[61,^, 74-7^ 
®’'Rb[61,^, 74, 7^ 
®’'Rb[^,_77J 
39K[^,22] 

®’'Rb[M, 78] 

^’^Rbf24. 25] 

Cs[79] 

®7Rb,[^,^ Cs[^ 

Cs[^ 


The early experimental studies of the Lieb-Linger gas 
with cold atoms were made in the laboratory by Bloch’s 
groupand Weiss’s group, [j^ In particular, the obser¬ 
vation of the ground state energy of the Tonks-Girardeau 
gas provides deep insights into understanding fermion¬ 
ization effect induced by a strong repulsive interaction, 
see Fig. , Loading the ®^Rb ultracold atoms into a 2D 
array of ID tubes, where the atoms were kept in the low¬ 
est energy state in the two transverse directions. Thus, 
the systems were realized in quasi-ID systems within ax¬ 
ial harmonic traps. The essential feature of the Tonks- 
Girardeau gas was observed through the ground-state en- 
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ergy Tid of such waveguided ensembles. 



Fig. 11. The ID ground-state energy Tid versus transverse 
confinement depth of the lattice. For the confinement po¬ 
tential Uq > 0 (or the effective interaction 7 1), the en¬ 

ergy Tid well presents the ground-state energy of the Lieb— 
Liniger gas in strong coupling regime within the local density 
approximation. [T^ 


The experimental measurement of the metastable 
highly excited state — the super Tonks-Girardeau gas 
was achieved by Haller in 2009. They made a new 
experimental breakthrough with a stable highly excited 
gas-like phase in the strongly attractive regime of bosonic 
Cesium atoms across a confinement-induced resonance, 
see Fig. □ This particular state was first predicted the¬ 
oretically by Astrakharchik et aZ.js^ Using the Monte 
Carlo method and by ANU group from the integrable 
interacting Bose gas with attractive interactions. ^ 
This model has improved our understanding of quantum 
statistics and dynamical interaction effect in many-body 
physics. It turns out that a highly excited state of gas¬ 
like gas could be stable as the interaction is switched 
from strongly repulsive into strongly attractive interac¬ 
tions due to the the existence of Fermi-like pressure. 32- 
This phenomenon has triggered much attention in 
theory. 



Fig. 12. The ratio of compress mode over the trapping fre¬ 
quency R = I ojp versus the interaction parameter . The 
squares show the experimental measurements in strong attrac¬ 
tive regime. The circles show the experimental data ranging 
from weak coupling to strong TG regimes. The black solid line 
stands for the exact result from the Lieb-Linger gas with a re¬ 
pulsion. The dashed lines present the theoretical data from 
the result. [ 30 ! Further study of the Tonks—Girardeau gas can 
be found in Ref. . 



Fig. 13. The local pair correlation function versus the effec¬ 
tive coupling constant. The solid line is obtained from the 
exactly solved model of Lieb-Linger gas at zero temperature. 
The symbols show the experimental data, [s^ 


In fact, many experiments have successfully demon¬ 
strated the confinements of atoms to one dimension 
by optical waveguides. Another particularly interesting 
example involves the measurement of photoassociation 
rates in one-dimensional Bose gases of ®^Rb atoms to de¬ 
termine the local pair correlation function (0) over 
a range of interaction strengths, see Fig. □ This ex¬ 
periment provides a direct observation of the fermion- 
ization of bosons with increasing interaction stren gth . 


It sheds light on the phase coherence behavior. [29|, 134 


42 I . l43l l81| At zero temperature, the local pair correla¬ 
tion is ~ 1 for the weakly interacting Bose gas 

and —>■ 0 as the system enters into the Tonks- 

Girardeau regime. 

As discussed in previous sections, the finite- 
temperature problem of the Lieb-Liniger Bose gas was 
solved by Yang and Yang in 1969.It turns out that the 
Yang-Yang thermodynamic equation is an elegant way to 
analytically access the thermodynamics, quantum Auc¬ 
tions, and quantum criticality. The Yang-Yang thermo¬ 
dynamics have been confirmed in the recent experiments 
through various thermodynamical properties and 

quantum fluctuations. 65l-l67j| A typical example is the 
measurement of the Yang-Yang thermodynamics on the 
atom chip, see Fig. „ 
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-400 0 400 -400 0 400 

X (^m) X (|im) 

Fig. 14. The in situ axial density profiles for the weakly 
interaction Bose gas of ®^Rb atoms at different temperatures. 
The solid lines show the result obtained from the Yang—Yang 
equations. The values of the chemical potentials are indicated 
by the harmonically trapping potentials. The experimental 
data show a good agreement with the theoretical prediction 
from the Yang—Yang equation. ISSU 

Moreover, recent experimental simulations with ultra¬ 
cold atoms provide promising opportunities to test quan¬ 
tum dynamics of many-body systems. In particular, 
nonequilibrium evolution of an isolated system involves 
transport and quench dynamics beyond the usual ther¬ 
mal Gibbs mechanism where the ground state and low 
lying excitati on p lay an essential role. The experimen¬ 
tal studyl?^ l7l| of thermalization of ID ensemble of 
cold atoms has led to significant developments in this 
field. [H, [s^ In these experiments, it was demonstrated 
that quenching the dynamics into the isolated systems 
can lead to non-thermal distributions if conserved laws 
exist. So far a generalized Gibbs ensemble is believed 
to present the non-thermal distributions in the isolated 
systems with conserved laws. The many-body density 
matrix is written as 



in terms of conserved quantities Xm- Here Z = 
Trexp ^rnXm^ IS the partition function. The La¬ 

grange multipliers Am acting for maximization of the en¬ 
tropy are determined by the associated conserved laws. 
FigurenShows the quantum Newton’s cradle that provides 
an insightful signature of such a generalized Gibbs ensem¬ 


ble. It shows that the ID systems with many-conserved 
laws do not approach the thermal equilibrium. 


Position (pm) 

-500 0 500 



0 0.5 1.0 


Normalized opticai thickness 

Fig. 15. The time series of absorption images of the first oscilla¬ 
tion cycle for initial average peak coupling strength ro = 1. The 
two groups of the cold atoms were confined in one dimension and 
initially separated by grating pulses. They evolved from time to 
time and collided twice in the centre of the harmonic ID trap in 
each full cycle. The oscillations in ID Bose gas last for a long time, 
even without approaching equilibrium. [70ll 

Recently, Langen et aZ.[7^ showed that a degenerate 
ID Bose gas relaxes to a state that can be described 
by such a generalized Gibbs ensemble. By splitting 
a ID Bose gas into two halves, they prepared a non¬ 
equilibrium system of ®^Rb atoms trapped in an atomic 
chip. They measured the local relative phase profile (p(z) 
between the two halves. It was shown that most of the 
experimentally reachable initial states evolve in time into 
the steady states which can be determined within a rea¬ 
sonable precision by far less than N Lagrange multipliers. 
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It was particularly interesting to see that the experimen¬ 
tal data of the reduced values can be well fitted with 
about 10 modes although there exist a much larger num¬ 
ber of conserved quantities in the system. This research 
further opens the study of the generalized Gibbs ensem¬ 
ble for the quantum systems out of equilibrium. 


OUTLOOK 


We have introduced a fundamental understanding of 
many-body phenomena in the Lieb-Liniger model. The 
exact results for various physical properties of the Lieb- 
Liniger model at T = 0 and at finite temperature were 
obtained by using the Bethe ansatz equations. In par¬ 
ticular, we have presented a precise understanding of the 
excitation modes, Luttinger liquid, quantum statistics, 
quantum criticality, correlations, and dynamics in the 
context of Bethe ansatz. In fact, there have been great 
developments in the study of the Lieb-Liniger model in 


Refs. [Ill 
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via various methods, such as field theory 


methods, m Luttinger liquid theory and bosoniza- 
tion, etc. 3a, l85j It was shown that the repulsive Lieb- 
Liniger Bose gas can be obtained as the nonrelativis- 
tic limit of the sinh-Gordon model. 8^-881 Moreover, 
the study of the non-thermal distributions for the iso¬ 
lated systems with many conserved laws has attracted 
much attention. In this scenario, the generalized Gibbs 
ensemble 83, U 9^ has been used to study the thermal- 
ization of the isolated systems. More recently, there has 
been growing interest on quench dynamics in terms of 
the generalized Gibbs Ensemble. |93M98| This research has 
been becoming a new frontier in cold atoms and con¬ 
densed matter physics. It turns out that the integrable 
systems of this kind thus provide a promising platform 
to advance the basic understanding of new quantum ef¬ 
fects in many-body physics, such as few-body problems, 
universal thermodynamics, universal contact, quench dy¬ 
namics, and correlation functions. These studies will 
further place mathematical theories of exactly solvable 
models into the laboratory for a wide range of physical 
phenomena. 
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